Signatures of Statistical Decay 
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ABSTRACT 

The partition of decay energy between the kinetic energy of reaction products and 
their Q-value of formation is obtained in a statistical derivation appropriate to highly 
C^l . excited nuclei, and is shown to be in a constant ratio. We measure the kinetic energy 

fraction, R = Y,Ekm/(^Eki n + EQo), over a wide range of excitation energy for well- 
defined systems formed in the 35 C1+ 12 C reaction at 35^4 MeV. Relationships between 
' excitation energy, charged-particle multiplicity, and intermediate-mass-fragment mul- 

. tiplicity, observed in this work and in recent experiments by a number of other groups, 

' follow from the derivation of the average kinetic energies and Q-values. 

in 
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q | 1. Introduction 

A number of scaling phenomena and correlations between observables have re- 
cently been discovered in the deexcitation of highly excited nuclei. These include: 

C ■ 

• the correlation of Njmf, the number of intermediate-mass fragments, with N c , 
^ ! the total number of charged products,© 

• the correlation of N IMF with Z bo ^ nd , the total amount of charge contained in 
fragments heavier than hydrogen,B 

• the approximate proportionality of Njmf to E t , the measured transverse energyi 

• the linear relationship between N c and T, the nuclear temperature^ and 

• the scaling of IMF multiplicity yield ratios with Hi 

Each of these observations has been taken as indicative of statistical decay. We 
here introduce a schematic derivation for the partitioning of the decay energy between 
the kinetic and the mass-excess degrees of freedom. Because of its relatively small 
width at moderate multiplicities, the ratio of the kinetic energy of the products to the 
total available decay energy is proposed as an event-by-event signature for statistical 
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Fig. 1: Emission energy (E k in + Qa) per nucleon in the moving frame plotted as a function of the 
velocity of the center of mass of the detected fragments for 35A MeV 35 C1+ 12 C. Points represent a 
massive-transfer simulation, filtered by the experimental acceptance; arrows indicate the velocities 
of the projectile and of the projectile-target center of mass. 



decay. Further consequences of the derivation include the dependence of the mean 
multiplicities, < N c > and < N IMF >, on excitation energy and the prediction of 
multiplicity yield ratios. The derived quantities are confronted with our experimental 
results for the decay of light nuclear systems with measured velocities, sizes, and 
excitation energies. 



2. Derivation of Statistical Observables 

The probability of emitting a particle with kinetic energy, Ekin, from a state of 
excitation energy, E, is 

P{E, Ekin,Qo) OC Ekin Ginv(Ekin) TFa ~~> (1) 

Pi{E) 

where pf and pi are the final- and initial-state level densities and <Jinv{E k i n ) is the 
cross section for the inverse reaction. The mean kinetic energy may be derived by 
integrating E kin with the probability distribution of equation (1): 



Jo Qo E kin P(E kin )dE, 
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In the limit of high excitation energy and negligible emission barriers, this gives the 
traditional! results for the mean kinetic energy, 



< E Un >= 2^/E/a = 2T, (3) 

and its variance, 
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2E/a = 2T 2 , (4) 



where a is the nuclear level-density parameter. These results are appropriate for 
neutron emission and for charged-particle emission from light systems, where Coulomb 
barriers are very low. The many exit channels accessible to states of high excitation 
energy permit one to approximate the density of available ground-state Q-values as 
a continuous function, f(Qo). A similar integral may then be performed for Qo, also 



giving a proportionality to yE/a, with the proportionality constant depending on 
/(Qo). Thus, 



< Qo >= 2^/E/a = 2T (5) 

if the density of available exit channels is a linear function of Qq. For systems in the 
mass range studied experimentally here, this is a reasonable approximation. Whatever 
the proportionality factor, the kinetic energy fraction, 

(EE kin + SQ ) ' ^ 

should then be a constant, independent of excitation energy, with the specific case of 
< Qo >= 2T, giving < R >= 0.50. 

For statistical processes, the observed width depends on the number of samplings 
of the parent distribution. In our case the number of samplings is the number, N, of 
detected charged particles. For < R >= 0.50, that width is 

°R = —?=■ (7) 
AV2N 

We now examine the relative rates for barrier-dominated and barrier-independent 
emission. If the probability per unit time for emitting n particles is approximately 
proportional to exp(—n < AE > /T), where < AE > is the average deexcitation 
energy per particle emitted, then the particle emission rate is 

/ dn n exp(—nAE/T) T 
< n >= — — ; . - — - oc 



/ dn exp(-nAE/T) AE' vy 

For light- ion emission from a light, highly excited nucleus, we neglect the Coulomb 
barrier and obtain 

< AE > = < E kin > + <Q >= 4T. (9) 
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Fig. 2: Energy sums and distribution widths Fig. 3: Centroids and widths of kinetic energy 

(bars) for 35,4 MeV 35 C1+ 12 C as a function of fraction, R, for 35 C1 + 12 C at 35A MeV. Data 

the number of detected charged particles. The for multiplicities 2 to 12 are plotted, 
solid line indicates the beam energy. 



However, for emission of heavier fragments, barrier effects may dominate, so that 
< AE >~ B. A fixed or generic value of B, previously demonstrated^ to be 
appropriate for IMF emission, would imply the rates for the two processes to be 
in the ratio, 

< Kimf > T/B 

oc — — — oc T. 10 

< n > T/4T v ; 

If the average deexcitation per particle is < AE >= 4T, the mean number of particles 
emitted would be < N >= E/4T, giving 



iVoc y/E, (11) 

and variance, 

o% = N/4. (12) 

From the ratio of rates, 

Njmf<x E. (13) 

If the proportionality of transverse energy to excitation energy is assumed, then equa- 
tions (11) and (13) consolidate the first four points of the introduction. 
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3. Experimental Determination of Source Properties 

A systematic comparison of the derived quantities with experiment requires the 
event-by-event determination of the size, velocity and excitation energy of the decay- 
ing system. One type of reaction which lends itself to the isolation and measurement 
of an excited source is the massive transfer mechanism!. A beam of 35 A MeV 35 CI 
ions, provided by Chalk River's TASCC facility, collided with a 2-mg/cm 2 carbon tar- 
get, producing a variety of reaction products. The kinematics of the reaction served 
to focus the reaction products into the 6° to 25° angular range of our detector array, 
a close-packed assembly of 40 phoswich countersi. All ions were identified by atomic 
number and isotopic distributions were measured with a set of high-resolution Si/Csl 
detector telescopes.0 Our selection of massive transfer events was facilitated by the 
thresholds and limited angular acceptance of the array, which reduced our sensitivity 
to target-like "spectator" matter. 

To ensure that no major component of an event went undetected, we analyzed 
only events in which at least 15 units of charge were identified. The velocity of 
the center of mass of all detected ions was computed, and the kinetic energy of 
each product was calculated within the moving frame. Based on the average mass 
excess deduced from the isotope distributions for each detected ion, T,Q was added 
to T,E kin for each event, and the resulting deexcitation energy per detected nucleon 
was plotted as a function of source velocity. Fig. 1 shows the relationship between 
emission energy per nucleon and source velocity, starting at projectile-like values of 
source velocity and zero emission energy in the lower right portion of the figure, and 
extending to compound- nucleus velocity and large emission energy in the upper left. 
The entire range of massive transfer phenomena, as previously observed at lower 
energies^, is evident in the figure. The points superimposed on the figure represent 
centroids of the distributions in energy and velocity for the simulated massive transfer 
of zero to twelve nucleons from the carbon target to the heavier projectile. Decay 
of the excited system is simulated!!! by a Monte Carlo event generator and filtered 
by the experimental acceptance. Quantitative agreement with the data has been 
demonstrated^ for simultaneous projection of the same set of calculated results on 
both the energy and velocity axes. 

A massive transfer reaction has the useful feature that the mass of the recoiling 
system may be deduced from its velocity, allowing a model-dependent estimate of 
total excitation energy. The validity of this estimate can then be tested by the 
consistency check demonstrated in Fig. 2, where the excitation energy, recoil energy, 
and deduced total energy are plotted as a function of the number of detected charged 
particles. Here, the multiplicity-dependence of the momentum transfer and energy 
deposition is evident. Significantly, the energy sums yield the initial projectile energy, 
indicating that the efficiency corrections have been properly applied. The behavior 
of any unobserved component must then be consistent with that of the detected 
fragments. 
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Fig. 4: Centroids and widths of multiplicity distributions for 35 CI + 12 C at 35^4 MeV, plotted as a 
function of reconstructed excitation energy. The proportionality constant for the solid line, k = 0.39, 
was chosen to fit the data. 



4. Comparison of Derived and Measured Statistical Observables 

Fig. 3 shows the kinetic energy fraction, R, plotted as a function of excitation 
energy for events in which two to twelve charged products were detected. The ex- 
citation energy attributed to each point is the mean for a given multiplicity. Filled 
circles are the centroids of the R distributions, with bars indicating the distribution 
widths, (Jr. The < R >= 0.50 value and associated widths from equation (7) are 
indicated by the dashed curves. To investigate the effects of a more comprehensive 
statistical treatment, a calculation was performed with the code, GEMINI 0, with 
average gamma-ray emission energies obtained from a modified version of PACE 0, 
and the results for three excitation energies plotted as open circles. The data are 
obviously in good agreement with both the schematic calculation and the statistical 
code results, though some bias due to detector acceptance may be apparent at the 
lowest multiplicities. The statistical nature of the decay is reflected in the R distribu- 
tions over the range of excitation energies from 2 to nearly 7 MeV per nucleon. The 
narrowness of the distribution supports the prediction of equation (7) and indicates 
that, even for relatively low multiplicities, R may be useful as an event-by-event indi- 
cator of statistical decay. Indeed, the technique has recently been applied to a series 
of projectile fragmentation experiments^, giving < R > values from 0.50 to 0.57. 

The mean number of charged fragments emitted should, according to equation 
(11), be proportional to the square root of the excitation energy. Fig. 4 shows 
the centroids and widths of the measured multiplicity distributions as a function of 
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Fig. 5: Natural logarithm of multiplicity yield ratios for 35 C1 + 12 C at 35^4 MeV, plotted as a 
function of E^ 1 / 2 . Heavy lines represent ratios based on "barrier-independent" calculations of N 
and c; light lines represent constant slope. 



excitation energy. The solid line represents a square-root proportionality, and the 
widths determined by equation (12) are indicated by the dashed lines. The centroids 
and widths predicted by our schematic statistical derivation are in good agreement 
with experiment for the higher excitation energies, but at lower energies, where the 
Coulomb barrier might impose upper limits on the number of charged particles, the 
experimental widths are narrower. 

Moretto et alB have demonstrated that for IMF emission, which may be dominated 
by barrier effects, the natural logarithm of the multiplicity yield ratio is a linear 
function of E~ 1 / 2 . We plot this ratio for emission of all ions in Fig. 5. Note that the 
data do not, except at the lowest excitation energies, have the linear dependence upon 
E~ l l 2 expected from the systematics of ref. 5. Based on the behavior of the widths 
in Fig. 4, it would instead be reasonable to expect agreement at higher excitations 
with multiplicity yield ratios predicted for barrier-independent emission. In this case, 
the assumption of gaussian multiplicity distributions gives 

1 Y(N) } kVE ' 1 ; 

where k is the proportionality constant of equation (11). These "statistical" ratios, 
indicated by the heavy lines in Fig. 5, do, in fact, approximate the data at high 
excitation energies. For lower excitation energies, the constant slopes of the barrier- 
dominated systematics may be more appropriate. 
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5. Conclusions 



• A schematic calculation showed, for E 3> B, that the partition of decay energy 
between kinetic energy of emission and Qo should be in a constant ratio for 
statistical decay. 

• The kinetic energy fraction, R = Y,E kin / (Y,E kin + HQ ), and its width were 
measured for a well-determined reaction mechanism. 

• The predicted mean and distribution width for R were observed in the data and 
reproduced in calculations with a well-known statistical decay code. 

• At sub-vaporization excitation energies, emission at the barrier has a different 
energy dependence than emission well above the barrier. 

• The relationship between the two processes is exemplified by many of the cor- 
relations observed between charged-particle production, IMF production, and 
excitation energy. 
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